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Hydrodynamic Approach to Vortex Lifetime in Trapped Bose Condensates 
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We study a vortex in a two-dimensional, harmonically trapped Bose-Einstein condensate at 
zero temperature. Through a variational calculation using a trial condensate wave function and 
a nonlinear Schrodinger Lagrangian, we obtain the effective potential experienced by a vortex at 
an arbitrary position in the condensate, and find that an off-center vortex will move in a circular 
trajectory around the trap center. We find the frequency of this precession to be smaller than the 
elementary excitation frequencies in the cloud. 

We also study the radiation of sound from a moving vortex in an infinite, uniform system, and 
discuss the validity of this as an approximation for the trapped case. Furthermore, we estimate the 
lifetime of a vortex due to imperfections in the trapping potential. 

PACS numbers: 03.75.Fi,03.65.Db,05.30.Jp,32.80.Pj 
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The prospect of creating quantized vortices in trapped Bose-Einstein condensed gases (EEC's) has been an intensely 
discussed and studied subject in the last few years m^- Despite the considerable interest in this area of BEC studies, 
many of the most fundamental questions have yet to be answered, such as those concerning the stability and lifetime 
of such a state P,H-n|. In this article, we study the properties of a vortex in a trapped BEC from a hydrodynamic 
point of view. We confine the discussion to two-dimensional systems at zero temperature, and furthermore employ 
the limit of a condensate which is large in comparison with the size of the vortex. 

We will thus only be concerned with a system whose properties can be described by a nonlinear Schrodinger 
equation. The system that we have in mind is the dilute Bose gas, which is governed by the Gross-Pitaevskii equation 
P,P| at temperatures sufficiently low that the system may be described as a superfluid. It should, however, be noted 
that this type of equation is applicable to a wider class of systems than just zero-temperature dilute gases ||lO| , |ll[] . 

The stability of a vortex in a BEC is limited by several factors. At finite temperatures, the vortex may be destroyed 
due to collisions with thermal excitations. [ |l2| This will not be the subject of this study. Second, the vortex may decay 
spontaneously even at zero temperature through the excitation of modes (or, equivalently, the emission of phonons) 
in the cloud y|7|]. Third, deviations from spherical symmetry in the trapping potential will also limit the lifetime of 
the vortex. The two latter processes will be the subject of this paper. 

The paper is organized as follows. Sections ||-|^ arc concerned with the motion of an off-center vortex and its decay 
through phonon emission. In section O we introduce the model and the trial assumptions for the density and velocity 
distributions, and analyze the motion of an off-center vortex. In Sec. Ill, the parallel case of a precessing vortex in 
an infinite, homogeneous system is analyzed, and the power loss due to the radiation of sound is calculated. This 
can be thought of as a "semiclassical" approximation to the trapped case considered here. In Sec. |^, the validity of 
this approximation is discussed and a lower bound for the lifetime of a vortex is arrived at. In Sec. M, we find the 
characteristic time for destruction of a vortex due to deviations from cylindrical symmetry in the trapping potential, 
and finally, in Section VI, the results are summarized and discussed. 



II. CIRCULAR MOTION 



In the case of a harmonic trapping potential, the equation for the condensate wave function i^if), whose squared 
modulus gives the superfluid density distribution p{r), reads 
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ijjir) = fiip{r), 



(1) 



where /i is the chemical potential and C/q — 47r?i a/m is the effective interaction potential, with a being the s-wave 
scattering length. Assuming the kinetic energy to be negligible, we obtain the so-called Thomas-Fermi approximation 
p3| for the wave function for a non-rotating cloud. 



/ ^2x1/2 

^TF{r^ = VP^[l--^j , (2) 

with an associated density distribution ptf = |'0tfP- Here, the central density po = ^^/Uo and the Thomas-Fermi 

1 /q 

radius R — (2ii/m(Jf) . For a two-dimensional system, the wave function ^[r) is normalized according to 

d^f\W)t^v, (3) 

where u is the number of particles per unit length. As a measure of the influence of the inter-particle interactions on 
the system's properties, we define the dimensionless parameter 

7 = va. 

The Thomas-Fermi approximation is a valid one when 7 is large, and in that limit we have for a two-dimensional 
system 

R = aosc27'/', (4) 

1 /2 

where Oosc — {fi/mujt) is the oscillator length. 

We now turn to the problem of a cloud containing a singly quantized vortex at the position tq. In the limit of large 
7, the density distribution of the cloud will not be appreciably affected by the presence of a vortex, except in a region 
whose size is comparable to the healing length ^, defined as 



h^ _ 1 

^~\2mpUo ~ (87rpa)i/2- ^ ^ 

The healing length gives the length scale over which the wave function for a vortex in a homogeneous Bose gas 
increases from zero to its bulk value |§,§,|4| . For an untrapped system, the density p is the value of the density far 
from the vortex core. In the case of a trapped system, p must be taken to be the local Thomas-Fermi density at the 
point r{) in the absence of a vortex, p(ro), thus defining a local healing length 

^(^0) = 7^ — ( \ M/9 = / ° , ^ (6) 
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where ^0 = ?(0) is the value of the healing length in the center. The velocity distribution in a Bose-condensed system 
is given by ti/m times the gradient of the phase of the wave function. For a positively oriented vortex in an infinite, 
uniform system, it is known to be 

m 
where (j) is the polar angle relative to the position of the vortex, which gives 

^ ,^, _ h z x (f~ fp) 
m \r ~ ?'o| 

z being the unit vector in the z direction. 

The velocity field is altered due to the boundary of the system and due to the spatially varying density. The 
existence of a boundary requires that the normal velocity vanishes there, and for homogeneous systems the recipe is 
to introduce a negatively oriented image vortex at the point fi = foR^ /rQ, giving the velocity field R 



'^The existence of sharp boundary is an artefact of the Thomas- Fermi approximation. Although this approximation never holds 
at the boundary, we do have that for sufficiently large clouds, p{R) is small, while Vp(7?) is nonnegligible, which justifies the 
neglect of the first term in the equation for stationary fiow, pV ■ v + v- V p = 0. Hence the radial velocity has to (approximately) 
vanish at r = R. 
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If the system has a density gradient, as in the present case, the condition for stationary flow, V • (pv) = 0, is not 
automaticaUy fulfilled. Writing the velocity field as 

w = Wo + u 1 , 

with vq taken from Eq. (M) , we get an equation for the correction vi : 

pV -vi+VQ-Vp + vi-Vp^ 0, (8) 

since the divergence of vq vanishes. An approximate solution to this equation, valid close to the center of the system, 
may be found by treating Vp as small, whereupon vi can also be expected to be a small correction, and the third 
term on the left-hand side of Eq. (ph can be discarded. We then have 

V • vi{r^ = /(r), 

where 
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with the boundary condition that the normal velocity vanish for r ~ R. Since we consider a Bose condensate, the 
velocity has to be a potential flow, vi{f) — V(/)i(r), and we have the equation for 0i, 

V2(/)i(r) = /(r) for r<i?, 

— = for r = it 

or 

whose solution is written in terms of the Green's function for the Neumann problem on a disk of radius R 

27r |r — 7^ 27r r* — \r\ 27r r 
as 

(/.i(f)=y'dV/(r')GAr(r',rO. (9) 

One can easily obtain higher-order velocity terms in Vp, if one writes 

V = Vi^ + Vi + V2 + ... 

One immediately finds that 



where w„ = V(/)„, n — 1,2, ... . When the density varies over a scale larger than the healing length, higher order 
corrections are small. We will be content in this paper to retain only the zeroth-order term. 

The energy per unit length of a two-dimensional system described by a nonlinear Schrodinger equation is 

£;[^,fo] ^j£r(^^ |V^|' + ^mcuy \i^f + ^ |^|^^ , (10) 

where we have made explicit the dependence of E on the vortex coordinate vq. The change in energy of the system 
due to the presence of a vortex only shows up in the kinetic-energy term, as long as its effect on the density profile 
is neglected. We shall denote this additional energy by UcS, since one may regard it as an effective potential for the 
vortex, depending on the vortex position vq. In terms of the velocity field v{f) it is written 

UMro)^^ fd'rpirHT^'- (11) 



Using the lowest-order approximation vq for the velocity, and employing the Thomas-Fermi wave function 
result is (cf. g) 
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To obtain this result, one needs to exclude the vortex core, of size ^(''o) around the vortex position, from the radial 
integral for the integrals to converge. The first term dominates for small rg and is independent of the details of the 
model. 

The above analysis shows that as long as the system is not subjected to an external rotational constraint, the vortex 
will experience an effective potential which decreases with the distance from the trap center, except in a small region 
close to the boundary, where [/off has a local minimum due to the unphysical behaviour of the Thomas-Fermi wave 
function at r = i?. This minimum is located a distance 6 from the edge, given by 

It is interesting to note that the healing length does not set the length scale here. Apart from constants of order 
unity, this "boundary thickness" is the same as the cut-off length S of Ref . ||l^ , which comes about when computing 
the kinetic energy of a Thomas-Fermi cloud. 

We now turn to the motion of the vortex in this simple model. We assume that the vortex coordinate tq may 
have a time dependence, but that all other parameters of the system remain stationary. Solving the time-dependent 
counterpart to the equation ffl) is equivalent to minimizing the action obtained from the Lagrangian |I3,n6| 



L['0, n),fo] =T[ip,fQ,fo]- E[-tJj,fo], 



where the kinetic term 



T[7/',fo,7^o] 



dV- 



.ih 



'dtp dip* 



A straightforward calculation, remembering that the gradient of the phase of ■0 is the velocity field, yields 

.2- (ro X fo) 



T: 



{v{ro) - v) 



(14) 
(15) 

(16) 



where vif^) = Svr /„" drrp(r) is the number of particles per unit length inside the circle of radius tq. The Euler- 
Lagrange equation for fo and fo will finally yield, for the radial (tq) and azimuthal ((/iq) components respectively, 

ro = 0; 

hdv/dr^ 

where we have defined F — —dUcs/drQ. We see that in this model, an off-center vortex executes an orbiting motion 
around the center with an angular frequency lu = c/jq- Since dv/drQ — 2nrQp{ro), we finally obtain 

2Trnrop{ro) 
In the case of a Thomas-Fermi profile, F is obtained by differentiating Eq. (0 
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where 
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which yields the final result for the frequency of precession of a vortex in a Thomas-Fermi cloud, 

h 
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9{ro/R). 



(17) 
(18) 

(19) 



The same expression has been obtained previously by different approaches |17 
allows a smooth generalization to incorporate the phonon effect. 



181. Nevertheless, the above treatment 



III. PHONON RADIATION BY A VORTEX IN AN INFINITE SYSTEM 

We now turn to the problem of a vortex in an infinite system, exercising (e. g. under the influence of an external 
force) circular motion. 

It has previously been shown ||l9|,^ how any homogeneous superfluid described by a nonlinear Schrodinger-type 
energy functional is equivalent to (2+l)-dimensional electrodynamics, with vortices playing the role of charges and 
sound corresponding to electromagnetic radiation. For a fluid, which in the absence of vortices has the density po 
(note that the use of this symbol is not the same as in the preceding section) , with the local fluid velocity t/(r, t) and 
density p{r,t), and (possibly) containing a vortex at the position ro(i), moving at a velocity foit) — Vy(t), we define 
the "vortex charge" q^ — —Tiyj2'KpQ/'m, "vortex density" Pv{r,t) = 5^'^\r — ?\)(i)), and the corresponding "vortex 
current" Jt,(r, i) — qyPv{'r,t)vy{r,t). The speed of sound is c = ^/Tjopofm. We then have the analogous Maxwell 
equations: 

V-5 = 0, 

V • e = 2-Kqyp^, 

Vxe+-— = 0, 
c at 

^ -. Ide 27r ^ 

V X 6 — - = — q^]^, 

c at c 

where we have defined 



^ 1271771 , ^ , ^ ^/ ^ X 

e{r,t) = \ p{r,t)z x v{r,t), 

Pa 



T< rhrm . ,^ , 

K'T.t) = W czp(r,t). 

V Pa 

The "no magnetic monopole" law is clear from the definition of b; the Coulomb law states how the presence of vortices 
create a rotational current, the Faraday law is equivalent to the continuity equation for the fluid, and the counterpart 
to Ampere's law derives from the Josephson-Anderson relation implied by the Euler equation. The energy of the 
system is 

E=^ /"dV(e^(f,t) + ?(f,t) 

and correspondingly the Poynting vector 

a(r,t) = ^e(r,i)x&(r,i). (20) 

Electromagnetic potentials, d{f, t) and (p(r, i), are defined in the usual way, and within a Lorentz gauge we recover 
the usual wave equations 

V'-i^')(^(r,i) = -27rp,(f,t), 
^2 1 5^\ ^.^ X 27r^,^ , 



which have the exact solution in (2+1) dimensions [pT| 

f{r,t) = - [dt'd^7^^^ ^ ' Pvi^,t'), 



t-t' - £^ 



^t-t'Y-\s^ 



e t-t 



I \r—r 



{r,t) = -- fdt'd^T^ / ; Jj ^i^^i')- (21) 



The step functions theta in the numerator is a feature pecuHar to two dimensions; the corresponding three-dimensional 
expressions contain a delta function. 

We wish to use the above formulation of vortex dynamics in order to find the energy dissipated from a circularly 
moving vortex due to the radiation of sound waves. Our aim is therefore to find the value of the Poynting vector a 
at large distances from the vortex. 

We first specialize to the case of one point particle at the position fo ■ Eq. (|2l|) becomes 



ip{f,t) = -qy / dt' 



it-t'- 4^ 

2 _ X2" 



^Jit-t') 



C J 1(4. ^i\2 _ X^ 



^{t-t'Y 



(22) 



where X = \f — ro(i')l- We now seek the solutions to these integrals in the limit of large r. To zeroth order in r^/r, 
X = r, which is independent of t' . For a charge exercising circular motion with frequency w at a radius ro, the velocity 
is 

Vy{t ) = v(~xsinujt +ycosLut). 

The integrals for the vector potential's x and y components can now be done exactly, yielding Bessel functions: 

rt — r/c 1,1 ' if 

' ' at SUV Lot 



ax[r,t) = 

c 



ay{r,t) 



{t-t'r-r± 

QyV /■*"'■/'= dt' COS LOt' 



qVTT 

"27 



No f — J amujt + Ja f — J 



cos Lot 



^{t-t'f 



qVTT 

"27 



-iVo — cosLot + Jo — I smwi 



\ c / 
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and finally 



a(r,t) = — 



Vv{t)No i^—j +ZX Vy{t)Jo i^—j 



We immediately obtain the "magnetic field" : 

J- TTqyLO 



2c2 



Vy{t)Ni (—) +ZX Vy{t)Ji {—\ 



V c 



\ c / 



Finally, we utilize the asymptotic formulas for the Bessel functions at large arguments, which yields 



(23) 



K"^ = Qv 



2ch 



Vy{t) X rsm — I + (z X Vy{t)) x rcos — 



(24) 



It is not necessary to calculate the electric field e{r) in order to obtain the Poynting vector, since at large r, the 
field is locally that of a plane wave, in which case we have 



a = -— 6 r. 

27r 



(25) 



On integrating around the circle with radius r we get the power radiated by a vortex exercising circular motion in an 
infinite system: 



rdOr ■ a 



■nqloov"^ 
4c2 



4c2 ■ 



(26) 



IV. LIFETIME OF A VORTEX 



In Sec. y we found that an ofF-center vortex performs a circular motion. In the preceding section we saw how, in an 
infinite and homogeneous system, such motion excites sound waves, which carry away energy from the vortex. In a 
trapped cfoud, the effect of such radiation would be that the vortex move outward towards regions of lower potential 
energy UcS^ro), until it finally escapes from the cloud [^,0. 



The application of the results of Sec. Ill on the trapped case may be thought of as a semiclassical approximation. 
One condition for this approximation to hold is that the precession frequency of the vortex match the attainable 
excitation frequencies of the cloud, as seen in Eq. (124), where the moving vortex excites sound waves which have the 
same frequency as the precession. 

This requirement, however, is not met in the present case. In a harmonically trapped cloud containing a vortex, all 
but one of the mode frequencies are greater than or equal to the trap frequency ujt plp2| ; the single low- lying mode 
is identical to an off-center displacement of the vortex |j,^,^ . 

Comparing the precession frequency with the trap frequency, we find, using Eq. (}^, 

^ _ 9{ro/R) 



87V2 (1 - ^) 



which is always less than one, except very close to the boundary. 

We conclude that the semiclassical approximation is never valid for a vortex in a harmonically trapped cloudy. This 
does not, however, necessarily imply that the vortex is stable; only that its lifetime is longer than that implied by the 
semiclassical approximation. 

The results of the two preceding sections can therefore be utilized to calculate a lower bound for the vortex lifetime. 
The power dissipated from the vortex by phonon emission, Eq. (Eq), is to be set equal to the rate of motion downhill 
the potential gradient: 

p^dE ^ Fdrp 
dt dt 

We note that P and F are functions of tq. Rearranging terms, we obtain the time Tp for the vortex to move from 
To = ^ to To = i? — (5: 

This quantity is a lower bound for the lifetime of a vortex originally residing at a distance S, from the trap center. 
The upper cutoff, i? — (5, is needed in order to avoid unphysical boundary effects, as discussed in connection with Eq. 
(O). The subscri pt p is introduced to indicate that this time scale is associated with the radiation of phonons. 
Inserting Eqs. (|l9|, <^ and (||) into @ we get 

_ IQm'^poUoR'^ 



where the dimensionless integral / equals 



l-S/R ,-^ 

dx- 



i/R a;(g(.T)) 



2 



with g{x) given by Eq. (Ilq). An exact result for the integral I is easily obtained numerically; the result is shown 
in Fig. 1, and will be discussed shortly. An estimate can be obtained by noting that the function g{x) for strong 
coupling is approximately equal to 2 ln(i?/^) -I- -^ over a large range of values of x, and that the lowest-order term in 
X dominates the numerator, whereupon one gets / « ln(i?/^)/(21n(i?/^) -I- i)^ and 

^ _ 16mVo{/oi?Mn(i?/C) 



^This is, in fact, the case for a larger class of trapping potentials, including all power-law potentials and the square well. 



Finally, we insert the Thomas- Fermi results (||), to cast the above result in terms of the parameter 7: 

1 128 7^/2 ln(4^) 



ut IT (ln(4V7) + i) 



1^2- 



(29) 
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FIG. 1. Lower-bound estimates for the lifetime of the vortex, multiplied by the trap frequency tut- The solid line shows the 
time scale Tp connected with phonon radiation, integrated numerically, while the dashed line is the quantity tm connected with 
broken rotational symmetry, with the parameter e set to 0.001. 



V. BROKEN ROTATIONAL SYMMETRY: CHANGE OF ANGULAR MOMENTUM 



The other factor which at zero temperature may limit the lifetime of a vortex is deviations from rotational symmetry 
of the trapping potential. Even very small irregularities in the magnetic or electric fields used to trap the condensed 
gases in experiments may affect the possibility to keep a vortex in the system. 

We therefore consider a BEC in a somewhat deformed trap. The total torque due to inhomogeneities in the density 
and the trap is 



/ 



M^ / d^fp{r)fx VV{f), 



(30) 



where p^f) and V{r) are the actual (not exactly cylindrically symmetric) density and external potential, respectively. 
A natural unit for measuring the torque would be the potential energy associated with the trapping potential, which 
has the same units: 



Uo: 



d'^fp{r)V{r), 



whose value for a harmonic-oscillator trap and a Thomas-Fermi density profile equals 

n 



Uo: 



12 



Pomu}^ R 



Writing 



M = M ■ z 



tU„. 



we can use e as a (approximately independent of coupling strength) measure of the relative distortion of the trap 
and the density profile. Assuming e to be approximately constant over time, the time scale for the destruction of an 
initially centrally placed vortex is 



where Lq is the value of the angular momentum for a system with a central vortex; this is equal to h times the number 
of particles v per unit length; and so the time for the vortex to move out is 

6h 

tm 



eniLufR^ 



VI. CONCLUSIONS 

We have arrived at two estimates for vortex lifetime: one connected with phonon radiation, and one due to broken 
rotational symmetry. The former, r^, is an increasing function of coupling strength 7, whereas the latter, tm, decreases 
with increasing 7. This leaves us with a "window" at moderate values of coupling strength, where both time scales 
are reasonably large. 

The lower-bound time scales Tp and tm are plotted as functions of coupling strength in Fig. 1. Both are to be 
considered as lower bounds, since the assumptions applied in deriving them are very pessimistic. The parameter e is 
taken to be 10~^, which is an upper bound for the trap deformations attainable experimentally p4| . Trap frequencies 
are often around 100 s^"'^ in experiments; thus Tp will be longer than one second (wich is a typical order of magnitude 
for condensate lifetimes) as long as 7 is greater than unity, and tm is longer than one second for all 7 < 5000. This 
leaves us with a large parameter range, easily attainable experimentally, for which a vortex at zero temperature can 
be considered as long-lived; considering the conservative assumptions made here, the actual region of vortex stability 
is probably much larger than this analysis indicates. 
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